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Гнтелектуальна система комп’ютерно! математики 
для високопродуктивних обчислень 


В статье рассматриваются И иллюстрируются на ряде примеров ИЗ различных классов задач вычислительной 
математики проблемы компьютерного моделирования при решении научно-технических задач. Решение 
проблем предлагается В рамках интеллектуальной системы компьютерной математики для высокопроизво- 
дительных вычислений: реализации компьютерной математики средствами интеллектуализации процесса 
исследования и решения задач. 

Ключевые слова: интеллектуальные системы, достоверность решений, параллельные 
вычисления, компьютерная математика, интеллектуальное численное программное 
обеспечение. 
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Введение 


Высокопроизводительная вычислительная техника является одним из основных 
средств научных и инженерных исследований. 

Вместе с ростом возможностей компьютеров для научных и инженерных исследо- 
ваний растут и проблемы их создания и эксплуатации. Увеличение числа процессоров в 
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параллельных компьютерах в этой ситуации будет означать существенное увеличение 
коммуникационных потерь и снижение их эффективности. Уже сейчас имеются су- 
щественные различия за счет коммуникационных потерь между максимальной и экс- 
плуатационной производительностями. Разработке компьютерных методов высокопро- 
изводительных вычислений (параллельных, распределенных, гибридных) посвящены 
работы [1-4]. 

Кроме того проблема достоверности компьютерных решений с ростом объемов 
решаемых задач на таких компьютерах также усугубляется. Известно, что в ряде 
случаев при решении научных и инженерных задач на компьютерах пользователи 
получают машинные решения, не содержащие физического смысла. Это происходит 
по многим причинам, но прежде всего из-за погрешности в исходных данных, раз- 
личия свойств математических и машинных моделей задач, различия арифметики и 
машинной арифметики и т.д. Проблема исследования достоверности компьютерных 
решений остается одной из практически важных. 

Другой, не менее важной, актуальной проблемой практической реализации 
высокопроизводительных вычислений является создание программного обеспечения 
уровня конечного пользователя — интеллектуальных программных средств, обеспе- 
чивающих общение с компьютером на языке предметной области и автоматизацию 
решения задачи на компьютере (алгоритмизация, программирование, выбор топологии, 
решение задачи в условиях приближенных исходных данных с анализом достовер- 
ности компьютерных решений). 

Решение перечисленных проблем предлагается в рамках концепции создания 
интеллектуальной системы компьютерной математики (ИСКМ): реализации компью- 
терной математики для высокопроизводительных вычислений средствами интел- 
лектуализации процесса исследования и решения задач. Ее применение позволит 
существенно перераспределить работы по постановке и решению задач между поль- 
зователем и компьютером по сравнению с традиционными технологиями, сократить 
сроки разработки приложений для решения научно-технических задач и повысить 
качество получаемых компьютерных решений, даст возможность решать новые 
научно-технические задачи и организовывать численные эксперименты, существен- 
но сокращающие средства и время разработок объектов современной техники. 


| Особенности компьютерной математики 


Математические модели, описывающие прикладные задачи, всегда содержат 
погрешности в исходных данных. Но в подавляющем большинстве случаев при ис- 
следовании математических уравнений предполагается неявно, что исходные данные 
задачи точны. Характерной особенностью математических моделей с приближен- 
ными исходными данными является то, что их математические свойства априори 
неизвестны. В пределах заданного уровня погрешности могут быть как совместные, 
так и несовместные задачи, как корректно, так и некорректно поставленные, как 
плохо, так и хорошо обусловленные. При этом машинная модель задачи, которую в 
конечном итоге и приходится решать на компьютере, всегда имеет приближенный 
по отношению к исходной задаче характер (из-за наследственной погрешности ис- 
ходных данных, из-за погрешности дискретизации, из-за погрешности получения 
(ввода) числовых данных о задаче в компьютер). 

Влияние погрешности исходных данных на точность полученных решений 
проиллюстрируем на простейших примерах. 
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Свободные члены систем уравнений 


100х, + 500х, =1700 — 100х, +500х, =1700 
15х, + 75,01х, =255 15х, +75,01х, = 255,03 


различаются в пятой значащей цифре. Решением первой системы является х, =17, 
х, = 0, авторой-— х, =2, х, =3. Отметим, что определитель обеих систем равен единице. 
Не всегда близость элементов матриц обеспечивает близость их собственных 


1 10* 
значений. Так, собственными значениями матрицы являются Я, =1, 4, =|, а 
1 
1 10 
матрицы 7 — ДА =2, А, =0. В этом случае из-за возмущения элементов 
0 1 


меняется каноническая форма матрицы от Жордановой клетки до матрицы простой 
структуры. 
А № . . 
Уравнение х’? —2х/° +1=0 имеет один действительный корень х!=1, а урав- 


й 


1 1 
нение хз —2х/6 +0,99 =0, отличающееся от предыдущего уравнения лишь погреш- 


ностью свободного члена = =10” , имеет два действительных корня х! = 1,771561, х= 
= 0,531441. 

Иллюстрацию таких примеров можно продолжить также из других классов за- 
дач вычислительной математики. 

Следующий пример показывает, что увеличение точности задания исходных 
данных далеко не всегда может приводить к уточнению результатов решения. Чтобы 
проиллюстрировать этот факт для примера рассмотрим матрицы 


4 0 = ° ^0 
Н = иН= Е 
оо 0 10 
Очевидно, что соответствующие псевдообратные матрицы 
1 — 1 
Нов рй у и НТ! = й у 
оо 0 10° 


будут иметь сильно отличающиеся элементы. 
Используя представление 


т 
+ _ Не Т 
И Е 
1=1 
где [; — собственный вектор матрицы Н, Я;— собственные значения, 


+= РИ. если д *0 
| 0, если А=0 | 


легко понять этот факт, поскольку функция А,’ имеет разрыв в точке А =0. Причем, 


как следует из представления Н’ ,‚ различие между Н*и Н ' будет тем больше, чем 
меньше погрешность в исходных данных. Выше приведены примеры плохообуслов- 
ленных и некорректных задач в условиях приближенных исходных данных. 

В этом случае чрезвычайно усложняет ситуацию тот факт, что большое ма- 
тематическое различие между матрицами полного и неполного ранга существует 
только в математически идеальном мире вещественных чисел. Поскольку действия 
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над матрицами проводятся с округлением, то это различие становится неопределен- 
ным. Таким образом, некоторая невырожденная матрица может стать в компьютере 
вырожденной. С другой стороны, очень вероятно, что вырожденная в действитель- 
ности матрица за счет погрешностей округлений будет превращена в близкую, но 
невырожденную. 

Решение проблемы состоит в том, чтобы в машинной среде определить свойства 
решаемой задачи и сформировать машинный алгоритм получения приближенного 
решения математической задачи как для корректных задач, так и для некорректных, 
как плохо, так и хорошо обусловленных. 

Анализ особенностей реализации компьютерной арифметики показал, что: 

— континуум всех вещественных чисел в компьютере аппроксимируется ко- 
нечным множеством конечных дробей (уже при вводе численных данных возникают 
ошибки округления); 

— феномен «машинного нуля» порождает ряд трудностей при реализации вы- 
числительных алгоритмов (любой современный компьютер имеет наименьшее поло- 
жительное число, которое может быть в нем представлено, и все числа, меньше по 
абсолютной величине этого числа, заменяются нулем); 

— арифметические операции на компьютере отличаются от математических: 
законы ассоциативности и дистрибутивности не выполняются ни на одном совре- 
менном компьютере, а законы коммутативности в операциях с плавающей запятой 
выполняются только при правильной процедуре округления. 

Таким образом, аксиоматика математики, в том числе вычислительной мате- 
матики, отличается от аксиоматики компьютерной математики. 

В ходе численного решения уравнений машинных моделей задач алгоритмами 
прямых методов, которые представлены в компьютерах в машинных кодах, в зави- 
симости от длины мантиссы машинного слова происходит накопление погрешностей 
вычислений. Влияние вычислительной погрешности на решение систем линейных 
алгебраических уравнений видно на примере решения системы алгоритмами методов 
Банча и Гаусса. 

Так, после ввода системы линейных алгебраических уравнений 

Ах =Б, 

где А и В имеют следующий вид: 

0.1348531574394 0.1878970588235 0.1909117647058 0.17792647058823 


0.1878970588232 0.262 0.265 0.247 
— 0.1909117647058 0.265 0.281 0.266 
0.1779264705882 0.247 0.266 0.255 
Ь=0.3516, 0.4887, 5105, 0.4818, 
а решение: 


х = 0.666216...е12, -0.401689...е12, -0.166554...е12, 0.979729...е11, 
в компьютере возникает система 

Ах = В, 
с матрицей А! и правой частью Ри, имеющих вид: 


0.1348531574394 0.1878970588235 0.19091176470588 0.17792647058823 
0.1878970588235 0.2620000000000 0.265 0.24699999999999 
0.1909117647058 0.2650000000000 0.28100000000000 0.26600000000000 
0.1779264705882 0.2469999999999 0.26600000000000 0.25500000000000 
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Ь: = 0.3516, 0.4887, 5105, 0.4818, 
точное решение которой: 
х! = 3.547...е12, -2.138...е12, -8.867...е12, 5.216...е11. 

Машинные решения этой системы как алгоритмом метода Банча 
Хванск = 2.810...е12, -1.694...е12,...-7.027...е11, 4.133...е11, 
так и алгоритмом метода Гаусса 
Хбаиз; = 3.164...е12, -1.908...е12,...-7.911...е11, 4.653...е11 
весьма далеки как от решения машинной модели задачи, так и от математического ее 
решения. 

Кроме того, при разработке вычислительных схем и программ возникает про- 
блема неоднозначной машинной реализации математических операций. 

Рассмотрим уравнение: 


+=) 
1(х)=0,001-100е` *’=0, 0<х<10, 
имеющее математическое решение 
х=5(6ш10—10)= 7,567. 
Для решения этого уравнения можно использовать метод простой итерации, 
реализуемый по формуле: 


хи =х, = 0,001+ и’ | 


так как выполнены достаточные условия сходимости этого метода 
р -®=) 
Ф'(х)=1-20е р 05-10: 


Решение, полученное на компьютере с ДЛИНОЙ машинного слова в четыре 


десятичных знака и с условием окончания итераций ыы] не | И пе ДЛЯ 
|х 1+= 
= =10° было получено за 3145 итераций: 
х=5,005, 

в то время как математическое решение есть 7,565, т.е. относительная погрешность 
машинного решения 33,8%. Увеличение длины машинного слова до 8 десятичных 
знаков по той же программе и с тем же условием окончания итераций для = =10* 
позволило получить решение за 33 237 итераций: 


х= 7,559 6341 


с относительной погрешностью 0,53%. 

Исходя из вышеизложенного, приходим к необходимости исследования мате- 
матических свойств компьютерных моделей задач. 

Например, для систем линейных алгебраических уравнений с целью выбора 
необходимого алгоритма решения и получения достоверных результатов целесо- 
образно исследовать: 

— корректность постановки математических задач, 

— определить или оценить число обусловленности сформулированной системы 
линейных алгебраических уравнений, 

— оценить погрешность машинного решения. 

Известно, что определяющим фактором при исследовании корректности систем 
линейных алгебраических уравнений является ранг матрицы и, в частности, для 
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квадратных матриц — исследование матрицы на сингулярность. Но использование 
классического определения ранга матрицы сводится, как известно, к вычислению 
определителя. Однако точно вычислить определитель матрицы на компьютере не- 
возможно из-за огромного объема операций вычислений и погрешностей машинной 
реализации, которые исказят истинный ранг матрицы. Сравнение вычисленного 
определителя с машинным нулем, очевидно, не даст правильного результата. Поэтому 
разработаны алгоритмы определения ранга матрицы в условиях приближенных ис- 
ходных данных, в том числе, и в условиях машинной погрешности. 

И, таким образом, для каждого класса математических задач с приближенными 
исходными данными возникает необходимость создания компьютерного инструмен- 
тария для исследования математических свойств машинных моделей задач, построения 
алгоритма их решения с учетом структуры и архитектуры компьютеров и оценки до- 
стоверности полученных результатов. В [1] представлены компьютерные алгоритмы 
высокопроизводительных вычислений. Вопросам достоверности получения компью- 
терных решений задач различных классов посвящены работы [5-10]. 


2 Интеллектуальное численное программное обеспечение 


Исследования в области компьютерной математики являются теоретической 
основой создания интеллектуального численного и прикладного программного обеспе- 
чения [11-17]. 

Под интеллектуальным программным обеспечением решения класса научно- 
технических задач будем понимать комплекс программ, позволяющий на языке 
предметной области сформулировать в компьютере задачу, автоматически иссле- 
довать свойства машинной модели задачи с приближенными исходными данными, в 
соответствии с выявленными свойствами и учетом математических и технических 
возможностей компьютера определить необходимое для решения задачи число про- 
цессоров и построить алгоритм решения, сформировать для решения задачи конфи- 
гурацию из процессоров параллельного компьютера, синтезировать программу парал- 
лельных вычислений, решить задачу, оценить достоверность полученного машинного 
решения и визуализировать результаты решения на языке предметной области. 

С функциональной точки зрения интеллектуальное программное обеспечение в 
автоматическом режиме реализует исследовательскую функцию и адаптивную на- 
стройку алгоритма, синтезированной программы и архитектуры компьютера на свойства 
решаемой задачи и получение компьютерного решения с оценкой достоверности. 

Алгоритмической основой интеллектуального программного обеспечения являются 
алгоритмы исследования и решения задач с приближенными исходными данными и 
оценкой погрешности получаемых компьютерных решений [2-4]. 

Возможности использования компьютеров для исследования математических 
свойств машинных моделей задач и их решения автоматически построенным ал- 
горитмом и программой были показаны в ходе выполнения проектов 1ЗРАК (Интеллек- 
туальное программное обеспечение для исследования и решения задач вычислительной 
математики на параллельных компьютерах) [18] и ГЗКОМ (Интеллектуальное про- 
граммное обеспечение для исследования и решения задач анализа прочности конст- 
рукций) [19], выполнявшихся для Немецкого центра по авиакосмическим полетам (ОГК). 

На рис. 1 представлен состав интеллектуального программного обеспечения 
прагзой, а также его взаимодействие с прикладным программным обеспечением. 
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Структурно Шшрагзой состоит из интеллектуального программного средства шрайоо| 
для автоматического исследования и решения задач и библиотеки интеллектуальных 
программ шрагИБ. 


Интеллектуальное Прикладное 
программное р 
обеспечение На 

Тпрагзойе обеспечение 
Тпрагкоо! Данные 
(интеллектуальное Решение отдельных 


программное средство Результаты 


решения задач 
вычислительной математики) 


прикладных задач 
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Прикладные 
программные средства 


Рисунок 1 — Состав шрагзой и его взаимодействие 
с прикладным программным обеспечением 


Интеллектуальное программное средство шрамоо| предназначено для иссле- 
дования и решения основных классов задач вычислительной математики, а именно: 

— систем линейных алгебраических уравнений; 

— алгебраической проблемы собственных значений; 

— уравнений и систем нелинейных уравнений; 

— систем обыкновенных дифференциальных уравнений с начальными условиями 
и реализует следующие возможности: 

— исследование математических свойств и решение задач с приближенными 
исходными данными; 

— автоматический выбор эффективных алгоритмов в соответствии с выявлен- 
ными свойствами решаемых задач, математическими и техническими характеристи- 
ками МИМО-компьютера; 

— исследование достоверности компьютерных решений; 

— реализация принципа скрытого параллелизма. 

Реализация принципа скрытого параллелизма обеспечивает: 

— автоматическое управление распределением и перераспределением инфор- 
мации между процессами; исключение эффекта Гайдна; 

— автоматическое определение необходимого для эффективного решения задачи 
числа процессов; 

— автоматическое построение эффективной топологии межпроцессорных связей. 

шрагюо| состоит из четырех отдельных компонент — каждая для исследования 
и решения задач одного из четырех названных выше основных классов задач вы- 
числительной математики. 

Функционально шрагзой является, с одной стороны, инструментом для автома- 
тического исследования и решения перечисленных выше основных классов задач 
вычислительной математики, а с другой стороны, параллельные программы библио- 
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теки шрагШЬ служат материалом (как геизе-программы) для создания прикладного 
программного обеспечения решения научных и инженерных задач. С точки зрения 
пользователя шрабоо! — продукт конечного пользователя, а шрагШЬ — инструмент 


пользователя-разработчика прикладных программных систем. 


На рис. 2 представлена блок-схема клиент-серверной архитектуры Шшрагоо]. 
Клиентская часть состоит только из диалоговой системы, а в серверную часть входят 
системы, обеспечивающие локальный и удаленный доступ пользователей к шрагоо|, 
а также системы, с помощью которых на ММО-компьютере проводится исследо- 


вание и решение задач с приближенными данными. 
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Рисунок 2 — Блок-схема архитектуры шра оо] 


Далее приведен фрагмент протокола автоматического исследования и решения 


задач на параллельном МИМО-компьютере Инпарком . 
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Мебсвоч: 
— Сайз$ е11т1пае1оп м1ЕБ рагЕ1а1 р1\уо®1п9 

5: МВ 5 
|1 ТНЕ МАТВТХ 15$ МАСНТМЕ-ЗТМСОТАВ 1!!! 
Мопрек оЕ ргосез5ог$: 4 
МмМебсвоач: 


— э1паи1аг уа1ае Чесотро$1Е1оп 

ОЕ а чепега1 мабг1х 
вт: 
СОТОТТОМ МАЗ САПГСОТАТЕР 


= 


Е1гзе 4 сошропепез оЕЁ зо1ае1оп (уесбог 1) аге: 
—-3.7747582837255322е-010 1.0000000000000031е+000 
3.8857805861880479е-010 3.6489927986770073е-010 
Тре уесгког (5$) оЕ зо1аЕ1оп аге зассеззЕа11у звогеа 10 ЕБе Е11е 
гези1 . ое 
Еггог езе1та®1оп5$: 4.99145е-08 
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— езе1таЕ1оп оЕЁ сопа1Е1опа1 папфег: 7.49316е+07 
— пабг1х гапКк: 999 

Мипрег оЕЁ ргосеззогз: 12 

Из протокола видно, что в качестве пробного алгоритма исследования задачи с 
помощью ШшрайМоо| выбран алгоритм метода Холецкого как наиболее экономичный 
алгоритм для решения систем с симметричными матрицами. Однако в процессе иссле- 
дования оказалось, что матрица не является положительно определенной, и для даль- 
нейшего исследования с помощью шраНоо| выбран алгоритм метода Гаусса. В процессе 
исследования матрица СЛАУ оказалась вырожденной. Для такой СЛАУ с использованием 
сингулярного разложения матрицы получено нормальное псевдорешение компьютер- 
ной задачи. Задача решена с оценками достоверности решения. В процессе исследо- 
вания и решения задачи автоматически устанавливалось оптимальное количество про- 
цессов для каждого алгоритма, строилась эффективная топология, в соответствии с 
алгоритмом решения распределялись (перераспределялись) данные между процессами. 

На базе интеллектуального программного обеспечения создан программный 
комплекс для автоматизации прочностных расчетов уникальных сооружений в граждан- 
ском и промышленном строительстве (НИр://уууу\и.Нга.сот.па). На Инпарком проанали- 
зирована прочность конструкций двухкорпусного 80-ти этажного офисного центра в 
Москве [20]. Показано существенное сокращение времени расчетов с использова- 
нием компьютерных технологий параллельных вычислений. 

Совместно с Институтом электросварки имени. Е.О. Патона НАН Украины ре- 
шена задача нахождения кинетики температурных напряжений и деформаций (на- 
пряженно-пластическая задача) при сварке титановых балок пола в конструкциях 
самолетов [21]. Использование технологии высокопроизводительных вычислений 
дало существенную экономию времени проведения расчетов. 


В 


ы 


Заключение 


Возросшие возможности вычислительной техники (высокая производительность 
и значительные объемы запоминающих устройств) дают возможность решать новые 
научно-технические задачи и организовывать численные эксперименты, существенно 
сокращающие средства и время разработки объектов современной техники. Отметим, 
что организация и стоимость натурных экспериментов требует затрат времени и 
стоимости на два — три порядка больше, чем организации численного эксперимента. 
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Теоретическое исследование возникающих математических и дискретных моде- 
лей является необходимым, но не достаточным условием получения достоверных 
компьютерных решений, так как из-за погрешностей свойства машинных моделей 
всегда будут отличаться от свойств дискретных моделей. 

Поэтому, интеллектуальные системы компьютерной математики для высоко- 
производительных вычислений является эффективным инструментом для автомати- 
ческого исследования, решения и анализа получаемых результатов задач инженерии 
и науки, максимально освобождая из этих процессов конечного пользователя, су- 
щественно повышая его производительность. 
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[иеШэеп! будет о Сотршег Майетайсх5 
юг Нйей-Регтотгтапсе Сотрийпе 


Тре рарег 4еа1$ уу гезоуше оРасбаа] ргоет сопсегише ргасйса| парететайоп оЁ 
Ые|-регюгтапсе сотрийпе Бу теапз оЁ сгеайоп оР зоЙ\уаге аё Ше еп4 изег’5 1еуе| — 
шеШоепЕ зузетз оЁ Ше сошрщег таетайс$ ргоу14те Бой Фе соттитшсаноп зв 
сотрщег ш Фе {егил$ ое заБ]есе агеа 1апоцазе ап4 ащютаноп оЁе ргоетл’5 зо оп 
сотприег (а|еогииайоп, ргозтатлиате, а свотсе оР‘юроозу, е зоуше оЁргоет ут 
фе арргохипае ша] даёа ‘осеег зу апау7те ое гепа Шу оРсотарщег зоаНоп5). 
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